Fusion of neutron rich oxygen isotopes in the crust of accreting neutron stars 
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Fusion reactions in the crust of an accreting neutron star are an important source of heat, and 
the depth at which these reactions occur is important for determining the temperature profile of the 
star. Fusion reactions depend strongly on the nuclear charge Z. Nuclei with Z < 6 can fuse at low 
densities in a liquid ocean. However, nuclei with Z = 8 or 10 may not burn until higher densities 
where the crust is solid and electron capture has made the nuclei neutron rich. We calculate the 
S factor for fusion reactions of neutron rich nuclei including O + O and Ne + Ne. We use 
a simple barrier penetration model. The 5* factor could be further enhanced by dynamical effects 
involving the neutron rich skin. This possible enhancement in S should be studied in the laboratory 
with neutron rich radioactive beams. We model the structure of the crust with molecular dynamics 
simulations. We find that the crust of accreting neutron stars may contain micro-crystals or regions 
of phase separation. Nevertheless, the screening factors that we determine for the enhancement of 
the rate of thermonuclear reactions are insensitive to these features. Finally, we calculate the rate 
of thermonuclear + ^"^O fusion and find that '^^O should burn at densities near 10" g/cm^ 
The energy released from this and similar reactions may be important for the temperature profile 
of the star. 

PACS numbers: 97.60.Jd, 26.60.+C, 97.80.Jp, 26.50.+X 



I. INTRODUCTION 

Nuclei accreting onto a neutron star undergo a variety 
of reactions. First at low densities, conventional ther- 
monuclear fusion takes place, see for example [ij. Next 
as nuclei are buried to higher densities, the rising elec- 
tron Fermi energy induces a series of electron captures 
0. Finally at very high densities, nuclei can fuse via 
pycnonuclear reactions. These reactions are induced by 
the quantum zero point motion [3]. The energy released, 
and the densities at which reactions occur, are important 
for determining the temperature profile of neutron star 
crusts. 

Superbursts are very energetic X-ray bursts from ac- 
creting neutron stars that are thought to involve the un- 
stable thermonuclear burning of carbon 0,111. However, 
some simulations do not reproduce the conditions needed 
for carbon ignition because they have too low tempera- 
tures An additional heat source, from fusion or other 
reactions, could raise the temperature and allow carbon 
ignition at densities that reproduce observed burst fre- 
quencies. 

Recently the cooling of two neutron stars has been ob- 
served after extended outbursts 0, Q . These outbursts 
heated the crusts out of equilibrium and then the cooling 
time was measured as the crusts returned to equilibrium. 
The surface temperature of the neutron star in KS 1731- 
260 decreased with an exponential time scale of 325 ± 
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100 days while MXB 1659-29 has a time scale of 505 ± 
59 days [8|. These cooling times depend on the ther- 
mal conductivity of the crust and the initial temperature 
profile. Comparing these observations of relatively rapid 
cooling to calculations by Rutledge et al. [P| and Shternin 
et al. suggests that the crust has a high thermal con- 
ductivity. However, if the initial temperature profile of 
the crust is peaked near the surface, then this peak could 
quickly diffuse to the surface and lead to rapid cooling. 
Therefore, cooling time scales are also sensitive to the 
initial temperature profile, and this depends on heating 
from nuclear reactions at moderate densities in the crust. 

Gupta et al. have calculated heating from electron 
capture reactions in the outer crust While they find 
more heating than previous works, they still find no more 
than 0.4 MeV per nucleon total heating from all of the 
electron captures on any mass number A system. Haensel 
and Zdunik have calculated pycnonuclear fusion reactions 
at great densities in the inner crust [ll|. However, if re- 
actions occur deep in the inner crust, most of the heat 
may flow in to the core instead of out towards the sur- 
face. As a result, there may be a smaller impact on the 
temperature profile of the outer crust. 

A low crust thermal conductivity, for example from 
an amorphous solid, could help explain superburst igni- 
tion. This could better insulate the outer crust and al- 
low higher carbon ignition temperatures. However, a low 
thermal conductivity appears to be directly contradicted 
by the observed short crust cooling times. Furthermore, 
our molecular dynamics simulations in ref. and fur- 
ther results we present in Section llVl find a regular crystal 
structure, even when the system has a complex composi- 
tion with many impurities. We do not find an amorphous 
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phase. These results will be discussed further in a later 
publication. We conclude that the thermal conductivity 
of the crust is high. 

If the thermal conductivity is high, one may need addi- 
tional heat sources, at moderate densities, in order to ex- 
plain superburst ignition. Although Gupta et al. find ad- 
ditional heating from electron captures to excited nuclear 
states, simple nuclear structure properties may provide 
a natural limit to the total heatin g fr om electron cap- 
tures fl^. Haensel and Zdunik [Til. Il3l| consider heating 
from pycnonuclear reactions using a simple one compo- 
nent plasma model. They find that fusion reactions may 
not take place until relatively high densities above 10^^ 
g/cw?. However, their use of a one component plasma 
could be a significant limitation. Fusion reactions depend 
strongly on the nuclear charge Z . Therefore, the reaction 
rate may be highest for the rare impurities that have the 
lowest Z, instead of for nuclei of average charge. 

In this paper, we go beyond Haensel and Zdunik and 
consider a full mixture of complex composition instead 
of assuming one average charge and mass. We focus 
on thermonuclear and pycnonuclear reactions at densities 
around 10^^ g/cvo?. This is near the base of the outer 
crust. Heat released at this density could be important 
for superburst ignition and for crust cooling times. Nuclei 
at this density are expected to be neutron rich. Further- 
more, the other nearby ions strongly screen the Coulomb 
barrier and greatly enhance the rate of thermonuclear 
reactions. 

We begin by describing the initial composition. This 
includes neutron rich light nuclei such as ^^O and ^^Ne. 
We calculate cross sections and S factors for ^''O -I- '^^O 
and ^^Ne -I- ^^Ne fusion using a simple barrier penetra- 
tion model. Note that the dynamics of the neutron rich 
skins of these nuclei can enhance the cross section over 
that predicted by our simple barrier penetration model. 
This is a very interesting and open nuclear structure 
question, see for example [141] . 

Next, we use classical molecular dynamics simulations 
to determine the structure of the crust and screening 
factors for the enhancement of thermonuclear reactions. 
There are many previous calculations of screening fac- 
tors for the one component plasma [l5j and for binary 
ion mixtures, see for example (l6j . However, we are not 
aware of any previous calculations for a crystal of a com- 
plex multicomponent composition. Finally, we calculate 
reaction rates and conclude that ^^O is expected to fuse 
at densities near 10^^ g/cva? while ^^Ne should react at 
densities near 10^^ g/cm'^. Heat from these reactions may 
be important for determining the temperature profile of 
accreting neutron stars. 



II. CRUST COMPOSITION 

We now describe our model for the composition of the 
crust. This is the same as was used in previous work on 
chemical separation when the crust freezes [ITi] . Schatz 



TABLE I: Abundance Xi (by number) of chemical element Z 
and average mass number {A). 
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Abundance (a;^) 




8 


0.0301 


24 


10 


0.0116 


28.8 


12 


0.0023 


36 


14 


0.0023 


42 


15 


0.0023 


45 


20 


0.0046 


62 


22 


0.0810 


66.06 


24 


0.0718 


74 


26 


0.1019 


76 


27 


0.0023 


77 


28 


0.0764 


80 


30 


0.0856 


89.35 


32 


0.0116 


96 


33 


0.1250 


99 


34 


0.3866 


102.61 


36 


0.0023 


106 


47 


0.0023 


109 



et al. have calculated the rapid proton capture (rp) pro- 
cess of hydrogen burning on the surface of an accreting 
neutron star [1], see also [l3|. This produces a variety of 
nuclei up to mass A k, 100. Gupta et al. then calculate 
how the composition of this rp process ash evolves, be- 
cause of electron capture and light particle reactions, as 
the material is buried by further accretion. Their final 
composition, at a density of 2.16 x 10^^ g/cw? , has forty 
% of the ions with atomic number Z — 34, while an ad- 
ditional 10% have Z = 33. The remaining 50% have a 
range of lower Z from 8 to 32. In particular about 3% 
is ^"^O and 1% ^*Ne. This Gupta et al. composition is 
listed in Table [H In general, nuclei at this depth in the 
crust are expected to be neutron rich because of electron 
capture. 



III. FUSION CROSS SECTIONS AND S 
FACTORS 

There is a great deal of experimental information on 
low energy fusion cross sections for light stable nuclei 
such as [H, 113 and ^^O For these nuclei, 

barrier penetration models work well (2l| . However, re- 
cently Jiang et al. discuss fusion hindrance at extreme 
sub coulomb barrier energies [l^ . Much less information 
is available for the fusion of very neutron rich light nuclei. 
We use a simple barrier penetration model to calculate 
fusion cross sections for ^^O and ^^Ne. We start with the 
Sao Paulo double folding potential Vf{t) [23l |. 

VF{r) = j d?rid\2Pi{ri)p2{r2)Vo&{vi-Y^-T). (1) 

Here pi and p2 are the densities of the two nuclei and 
Vo = -456 MeV-fm^. Next tunneling through the 
Coulomb barrier is calculated in a WKB approximation 
including some nonlocality effects [1, [2^ . For simplicity 
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TABLE II: Astrophysical 5* factors for low energy fusion re- 
actions versus center of mass energy E. 



(MeV) (MeV-barn) (MeV-barn) 



1 


1.7 X 10^^ 




1.5 


1.1 X 10^^ 


1.0 X 10' 


2 


7.4 X lO^"* 


7.2 X 10' 


3 


3.2 X 10^"* 


3.4 X 10' 


4 


1.3 X lO^"* 


1.6 X 10' 


6 


1.9 X 10^3 


3.4 X 10' 


8 


2.5 X 10^^ 


6.4 X 10' 



we assume Wood Saxon densities with radius parame- 
ter R = 1.31^1/3 - 0.84 fm and difTuseness a = 0.58 fm 
[1, [2^ . These parameters reproduce the measured cross 
sections for ^^O+^^O fusion. Our results for the fusion 
cross section (t{E) at center of mass energy E are ex- 
pressed as the astrophysical 5 factor, S = Ea{E)e'^'^^ 
and collected in Table [TTl Here the Gamow penetration 
factor is ?7 = ZiZ2e^(/x/2£^)^/^, the nuclei have charges 
Zi and Z2, and ^ is the reduced mass. Our S factor 
for ^'*0-|-^''0 is over eight orders of magnitude larger 
than that for ^^0-1-^^0. We have also calculated 5* for 
^'*0-|-^'*0 using relativistic mean field densities calculated 
with the NL3 interaction [24*1 . This yields S that is only 
slightly higher than the calculation with Wood Saxon 
densities. In addition, Casques et al. [1^ have calcu- 
lated S for ^"'O-I-^^'O using both a fermionic molecular 
dynamics model and the Sao Paulo model of Eq. [1] and 
they obtain similar results. 

These barrier penetration results may provide lower 
limits for the cross sections. Dynamical effects, not in- 
cluded in Eq. [1] can increase the cross section. Indeed, 
low energy cross sections are observed to be larger than 
simple barrier penetration estimates for heavier stable 
nuclei p6j . The extended neutron skin of very neutron 
rich nuclei presents a very interesting special case for low 
energy fusion reactions. The dynamics of the easily po- 
larizable skin can increase the cross section. For exam- 
ple, a neutron rich neck could form between the nuclei 
decreasing the Coulomb barrier. The dynamics of the 
skin, and its effects on low energy fusion, are very impor- 
tant nuclear structure questions that should be studied 
further with radioactive beams. For example, it should 
be possible to measure low energy fusion cross sections 
for beams of a neutron rich O isotope colliding with a 
stable light target such as ^^O or ^^C. 



IV. MOLECULAR DYNAMICS SIMULATIONS 
OF CRUST STRUCTURE 



al. [231 presented a phenomenological formula for reac- 
tions in a one component plasma (OCP) that is valid for 
all regimes of density and temperature. This has been 
extended by Yakovlev et al. to multicomponent plasmas 
(MCP) pj. 

In order to calculate pycnonuclear reactions in a mul- 
ticomponent system one needs to understand its state. 
Monte Carlo simulations [1^ of the freezing of a classi- 
cal OCP indicate that it can freeze into imperfect body 
centered cubic (bcc) or face-centered cubic (fee) micro- 
crystals. Unfortunetly not much has been published on 
the freezing of MCP. In an earlier work flT\ , we calculated 
chemical separation upon freezing of our MCP composed 
of rp process ash. We found, based on large scale molecu- 
lar dynamics (MD) simulations, that chemical separation 
takes place. The liquid phase is greatly enriched in low 
Z elements compared to the solid phase. We found that 
the solid phase formed a regular lattice where the charge 
of a given lattice site was more or less random. However, 
we did not study the structure of the solid in detail, aside 
from its average composition. 

There are many possibilities for the state of a cold 
MCP [2^. It can be a regular MCP lattice; or microcrys- 
tals; or an amorphous, uniformly mixed structure; or a 
lattice of one phase with random admixture of other ions; 
or even an ensemble of phase separated domains. We per- 
form classical MD simulations to explore the state of our 
MCP solid. The electrons form a very degenerate rela- 
tivistic electron gas that slightly screens the interaction 
between ions. We assume the potential Vij{r) between 
the ith and jth ion is. 



V^j{r) 



(2) 



where r is the distance between ions and the electron 
screening length is Ae = 7r^/^/2e(37r^ne)^/^. Here ng is 
the electron density. 

To characterize our simulations , we define Coulomb 
coupling parameters Tj for ions of charge Zj, 



Z}e^ 



(3) 



with T the temperature and is the radius of a sphere 
containing Zj electrons (the ion sphere radius). 



\4TTnpJ 



(4) 



The average coupling parameter P for the MCP is. 



aT 



We now consider fusion reactions in a dense medium. 
There have been many calculations of the strong screen- 
ing enhancement of thermonuclear reactions and of py- 
cnonuclear reactions. For example, recently Casques et 



where Xj is the abundance (by number) of charge Zj 
and the overall ion sphere radius a is a = (3/4?™)^/'^ 
with n = Ue/jZ) the ion density. The OCP freezes at 
r = 175. In ref. Il7l| we found that the impurities in our 
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TABLE III: Simulation Parameters. The temperature T, 
Coulomb parameter F, total simulation time t, average charge 
(Z), and average Z^^^ for the two simulations. Each simula- 
tion is at a density n = 7.18 x 10" fm"^ (1 x 10^^ g/cm^) 
and involves 3456 ions. 

Run r(MeV) F t(fm/c) {Z) {Z''^^) 



rpcrust-01 0.325 
rpcrust-02 0.35 



261.6 2.4 X lO'' 29.3 285.8 
242.9 1.6 X lO'' 29.3 285.8 



MCP lowered the melting temperature until T = 247. 
Finally, we can measure time in our simulation in units 
of one over the plasma frequency ojp, 



„ ZH'Ke^Xtn\i/2 



(6) 



higher temperature T = 0.35 MeV to possibly speed the 
diffusion of O ions. The total simulation time is 1.6 x 10^ 
fm/c. Note, that this run is ongoing and results for longer 
simulation times will be reported in a later publication. 
Figure [2] shows the final configuration of the 3456 ions. 
Now the system involves a single body-centered cubic 
(bcc) crystal. The O ions are not uniformly distributed. 
Instead O is strongly enriched in local regions. This 
simulation has F = 243. This corresponds to a signifi- 
cantly lower temperature than the melting temperature 
of a pure OCP which corresponds to F w 175. However, 
the impurities were found in ref. [l7j to lower the melting 
temperature till F « 247. Therefore the O rich regions 
in rpcrust-02 may be related to the formation of a bulk 
liquid phase which was found in ref. fv?\ to be greatly 
enriched in O. 



where Mj is the average mass of ions with charge Zj . 

To explore possible states for the multicomponent 
plasma we perform two molecular dynamics simulations. 
The initial conditions of these simulations are similar to 
those in We start by freezing a very small system 
of 432 ions. Here the ions were started with random 
initial conditions at a high temperature T and T was re- 
duced in stages (by re-scaling velocities) until the system 
freezes. For the first simulation run, called rpcrust-01 in 
Table IHIl we place four copies of this 432 ion solid in a 
larger simulation volume along with four copies of a 432 
ion liquid configuration. This 3456 ion configuration is 
evolved at a lower temperature until the whole system 
freezes. Finally, we evolve the 3456 ion solid at a refer- 
ence high density of n = 7.18 x 10"^ fm"'^ (or 1 x lO^'^ 
g/cm'^) and a temperature of T = 0.325 MeV for a to- 
tal simulation time of 2.4 x 10^ fm/c (8.9 x 10^ uj~^). 
This density and temperature correspond to F = 261.6. 
Evolution was done using the velocity verlet algorithm 
[3^ using a time step of At = 25 fm/c for a total of 
9.6 X 10^ steps. This took about 2 months on a sin- 
gle special purpose MDGRAPE-2 ^] board. The sim- 
ulation results can be scaled to other densities n' and 
temperatures T' that also correspond to F = 261.6. IE 
n'^/^/T' = (7.18 X 10-5fm"^)i/V(0.325MeV). 

The final configuration for run rpcrust-01, see Table 
mil is shown in Fig. [1] after a simulation time of 2.4 x 
10^ fm/c. The system is seen to be composed of two 
micro-crystals of different orientations. This is similar 
to the micro-crystals found in ref. [1^ upon freezing 
a one component plasma. In Fig. [T] we highlight the 
positions of the ^^O ions (as small red spheres). These 
ions are located both in the crystal planes and in between 
them. The O ions are not spread uniformly throughout 
the volume but there is a tendency for them to cluster. 
This will be discussed in more detail below. 

Given the micro-crytals for run rpcrust-01, we per- 
formed a second simulation, labled rpcrust-02 in Table 
mil with different initial conditions. Here eight copies of 
a 432 ion solid configuration were placed in the larger 
simulation volume. The system was evolved at a slightly 



FIG. 1: (Color on line) Configuration of the 3456 ion mixture 
in run rpcrust-01 after a simulation time of 2.4 x 10^ fm/c at 
r — 261.6. The small red spheres show the positions of ^'^O 
ions, while ions of above average Z are shown as larger blue 
spheres. Finally, ions of below average Z (except for O) are 
shown as small white spheres. The upper and lower halves of 
the figure show two micro-crystals of different orientations. 



We now use these simulation results to calculate the ef- 
fective screening potential (r) provided by all of the 
other ions. This greatly enhances the rate of thermonu- 
clear fusion reactions of two charge Zi ions. 



Zfe' 



(7) 



The radial distribution function gair) gives the proba- 
bility to find another ion of charge a distance r away 
from a given charge Zi ion. This is normalized to one at 
large distances gair — s- 00) — > 1. We calculate ga from 
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enriched in low Z ions. 




FIG. 2: (Color on line) Configuration of the 3456 ion mixture 
in run rpcrust-02 after a simulation time of 1.6 x 10^ fm/c at 
r — 242.9. The small red spheres show the positions of ^''O 
ions, while ions of above average Z are shown as larger blue 
spheres. Finally, ions of below average Z (except for O) are 
shown as small white spheres. The concentration is seen 
to be enhanced in sub-regions. 



our simulation by histograming relative distances. Fig- 
ure [3 shows gu{r) for = 8 (O), 10 (Ne), 22 (Ti), and 
34 (Se). Note that because of the low Ne abundance, our 
gu results for Zi = 10 are based on the positions of only 
40 ions! Therefore we caution that our Nc results may 
have large finite size and or statistical errors. 

Selenium is the dominate species. Therefore the Se 
ion locations largely determine the bulk structure of the 
crystal lattice. The distance between peaks in gu for Se 
reflects the lattice spacing. The radial distribution func- 
tion for Ti closely follows that for Se at large distances. 
This shows that Ti, for the most part, occupies the same 
lattice sites as Se. However the first peak in gu for Ti oc- 
curs at smaller distances than the first peak for Se. This 
reflects the smaller ion sphere radius aj Eq. [5] for the 
lower charged Ti because the coulomb repulsion between 
two Ti ions is smaller than that between two Se ions. 

The radial distribution functions for O and Ne show 
large peaks at small distances. This may reflect a ten- 
dency to replace a single Se ion with a cluster of multiple 
low charge O or Ne ions. The radial distribution function 
for O does not have large dips between the peaks in gu for 
Se. This shows that O also occupies positions in between 
the lattice planes. Finally, gu for O and Ne is larger than 
one at intermediate and large distances. This shows that 
the the low Z ions are not uniformly distributed. Instead 
they preferentially cluster in sub regions that are greatly 



FIG. 3: (Color on line) Radial distribution functions gii{r) 
for ions of charges Zi = 8(0), lO(Ne), 22(Ti), and 34 (Se). 
Results from run rpcrust-01 at a simulation time of 2.4 x 10^ 
fm/c have been scaled to a density of 3.6 x 10^" g/cm'^ and 
T = 0.05 MeV. 



The radial distribution functions for run rpcrust-02 are 
compared to run rpcrust-01 in Fig. ID Results for Se are 
similar. However, gui^) for O, at intermediate and large 
distances, is even more enhanced for rpcrust-02 than for 
rpcrust-01. Presumably, the slightly higher temperature 
of rpcrust-02 enhances phase separation into regions that 
are enriched in low Z ions. 

We now use these gu results to calculate reaction rates. 
Strong ion screening enhances the rate of thermonuclear 
fusion by a factor F, see for example [28j . 

F = exp(/iJ(r = 0)) . (8) 

Here h^{r) = — Wg^j(r)/T. For simplicity, we neglect the 

dependence of F on h^{r) for r ^ 0, see [sO]- Unfortu- 
nately, it is hard to get good statistics on gu (r) for small 
r. Therefore, one must extrapolate our MD results for 
r « flj to smaller r. The form of h{r) for small r is known 
for the OCP. We generalize the OCP expression for h{r) 
in ref. ^31i] to the multi-component case and assume for 
r < l.Soj, 

h^r)^hl^lT,{^)\0.0277T,{^y. (9) 

We fit Eq. [9] to our MD results for guir) over the range 
where we find a nonzero guif) and r < 1.5aj and extract 
values for = h^{r = 0). The enhancement of the 
thermonuclear rate is then F = Exp(/iq). Our values for 
/ig are collected in Table ITVl These values are averages 
of five values of /ig calculated at t=2.0, 2.1, 2.2, 2.3, and 
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FIG. 4: (Color on line) Radial distribution functions gii{r) for 
ions of charges Zi = 8 (solid lines) and 34 (dashed lines) at a 
density of 3.6 x 10^° g/cm^. Heavy lines show run rpcrust-02 
results at a simulation time of 1.6 x 10^ fm/c while thin lines 
show run rpcrust-01 results at £ = 2.4 x 10^ fm/c. 



2.4 X 10^ fm/c for run rpcrust-01 and at t=1.2, 1.3, 1.4, 
1.5, and 1.6 x 10^ fm/c for run rpcrust-02. At each time, 
we calculate /iq from the average of 25000 configurations 
separated by 250 fm/c. A simple analytic formula for hg, 
based on the linear mixing rule for the MCP is, see for 
example [28|, 



h\r = 0) « 1.0754rj 



(10) 



We find good agreement between our MD simulation re- 
sults and Eq. pH)) even at our large value of F = 262. 
Furthermore this is true for both rpcrust-01 and rpcrust- 
02 runs. For the large Z ions Ti and Se, the agreement 
between our MD results and Eq. ([TO]) is better than 1%. 
For the lighter ions O and Ne, our MD results are slightly 
larger than Eq. pO)) . Again we caution that our MD re- 
sults for Ne are based on only 40 ions. Therefore we focus 
on O. Our enhancement for O may reflect the clustering 
of O into sub-regions as shown in Fig. [21 

In any case, the overall agreement between our MD 
results and Eq. (jTU]) is very good. It is a major re- 
sult of this paper that features in the crust, such as the 
micro-crystals in Fig. [1] or the phase separation in Fig. 
[21 do not appear to be important for the screening of 
thermonuclear reactions. Therefore, for the rest of this 
paper we will simply use Eq. (jlO[) to describe screening. 
Note that these features in the crust may be much more 
important for calculating pycno-nuclear reaction rates at 
very high densities. This will be explored in later work. 



TABLE IV: Screening potential h^r = 0) for O, Ne, Ti, and 
Se ions {Zj = 8, 10, 22, 34) from our MD simulations. Also 
listed is an analytic approximation h-' {r = 0) = 1.0754rj, 
Eq. (|10|l and the difference between our MD results and the 
analytic approximation. 

Run r O Ne Ti Se 

rpcrust-01 261.6 32.5 48.6 169.6 350.0 

analytic 261.6 31.5 45.7 170.0 351.2 

MD - analytic 261.6 1.0 2.9 -0.4 -1.2 

rpcrust-02 242.9 30.3 43.8 158.1 325.5 

analytic 242.9 29.2 42.4 157.8 326.1 

MD - analytic 242.9 1.1 1.4 0.3 -0.6 



V. REACTION RATES 

We now calculate the rate of thermonuclear ^^0+^*^0 
fusion including the effects of strong screening. The re- 
action rate per O ion R/n is given by the well known 
formula, see for example [28j . 



n \ 3/i 



L/2 S{E. 



T 



(11) 



Here r = (277r^^Z^^Q;^/2T)^/'^, ^ is the reduced mass 
and -Epfe = tT/Z. We consider a typical temperature 
of 5.8 x 10^ K [3]. The reaction rate versus density is 
collected in Tabic [Vl These results assume the S factor 
from Table [in the ion screening from Eq. (jlOp . and the 
number fraction of ^"^O is 0.1. Oxygen will burn at a 
density where the reaction rate per ion, times the time 
for a fluid element to be buried to a given density, is one. 
It can take of order 1000 years (3 X 10^° s), depending 
of course on the accretion rate, for a fluid element to be 
buried to these densities. We conclude from Table [Vl that 
^^O will burn at a density near 10^^ g/cm"^. 

The neutron rich ^^O has a lower fusion rate than ^^O 
because of its larger reduced mass. However, most of this 
reduction in rate is compensated for by a much larger S 
factor. Our S factor in Table [TT] is over eight orders of 
magnitude larger than the S factor for ^^O+^^O fusion. 
As a result the thermonuclear rate Eq. pT|) for ^''O is 
only slightly smaller than the rate for -'^^O fusion. Indeed, 
if there is a significant enhancement in 0-1-^*0 fusion 
because of the dynamics of the neutron rich skin, the ^''O 
rate will be larger than the ^^O rate. 

Table [Vl also lists the ratio of temperature to ion 
plasma frequency, Eq. Strictly speaking the ther- 

monuclear fusion rate, Eq. pip , is only valid for T > ojp. 
We see that the rates in Table [Vl involve an extrapola- 
tion of Eq. (fTTl) to T slightly below ojp. At these tem- 
peratures, there will be some quantum corrections to Eq. 
(jlip . However, the thermally enhanced pycnonuclear fu- 
sion rates in [2^ suggest that Eq. ([11]) is not wildly wrong 
at these temperatures. Therefore, quantum corrections 
should not change our conclusions that ^"^O will burn at 
a density near 10^^ g/cm^. 
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TABLE V: Reaction rate for ^*0+^''0 fusion versus density 
p. The coulomb parameter is F, while the ratio of the tem- 
perature to the ion plasma frequency is T/ojp, and the fusion 
rate per ion is Rjn. 

p (g/cm3) r T/iUp R/n (s"^) 

10^" 170.0 2.2 6.9 X 10"^^ 

4 X 10^" 269.9 1.1 4.5 x 10"" 

10" 366.3 0.68 1.2 x 10"" 

2 X 10" 461.5 0.49 2.4 x 10"'' 

We also consider ^^Ne+^^Ne fusion since ^^Ne is the 
next heavier nucleus in our rp process ash, after ^^O. 
However, because of the larger charge Z — 10, we find 
that ^*Ne will not burn until higher densities where the 
plasma frequency is much larger than the temperature. 
Therefore ^^Ne will not burn via thermonuclear fusion. 
Instead, it will burn via pycnonuclear or thermally en- 
hanced pycnonuclear fusion. Using the pycnonuclear 
rates in ref. [1^ along with our S factor from Table 
im we estimate that ^^Ne will burn at densities near 10^^ 
g/cvD? (for temperatures near 5.8 x lO^K). 

The fusion of ^^O releases 0.52 MeV per nucleon while 
^^Ne fusion releases 0.64 MeV per nucleon. These en- 
ergies are larger than the total heating from all of the 
electron captures considered by Gupta et al. [1]. There- 
fore, these and other related fusion reactions may be an 
important source of heat in the crust of accreting neutron 
stars. 



VI. SUMMARY AND CONCLUSIONS 

Fusion reactions in the crust of an accreting neutron 
star are an important source of heat, and the depth at 
which these reactions occur is important for determin- 
ing the temperature profile of the star. Fusion reactions 
depend strongly on the nuclear charge Z . Nuclei with 
Z < 6 can fuse at low densities in a liquid ocean. How- 
ever, nuclei with Z = 8 or 10 may not burn until higher 
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densities where the crust is solid and electron capture has 
made the nuclei neutron rich. 

In Section lHll we calculated the S factor for fusion reac- 
tions of neutron rich nuclei including ^'^O -f ^'^O and ^^Ne 
+ ^^Ne. We used a simple barrier penetration model. We 
find that S for ■^^O is over eight orders of magni- 

tude larger than that for ^^0-1-^^0. The S factor could 
be further enhanced by dynamical effects involving the 
neutron rich skin of ^^O. For example, the skins of the 
two nuclei could deform to form a neck that would re- 
duce the Coulomb barrier. This possible enhancement in 
S should be studied in the laboratory with neutron rich 
radioactive beams. 

In Section IIVI we modeled the structure of the crust 
with molecular dynamics simulations. We find that 
the crust of accreting neutron stars may contain micro- 
crystals or regions of phase separation. Nevertheless, the 
screening factors that we determined for the enhance- 
ment of the rate of thermonuclear reactions are insensi- 
tive to these features. Finally, we calculated in Section IVl 
the rate of thermonuclear ^^O + ^^O fusion and find that 
^^O should burn at densities near 10^^ g/cvo?. This is a 
lower density than some previous estimates. The 0.52 
MeV per nucleon energy released may be important for 
the temperature profile of the star. In future work, we 
will use our molecular dynamics results to study other 
properties of the crust such as its thermal conductivity. 
In addition, we will use these MD results to calculate py- 
cnonuclear reaction rates for the fusion of ^^Ne and other 
heavier nuclei. 
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